Planar QED at finite temperature and density: Hall 
conductivity, Berry's phases and minimal 
conductivity of graphene 



C.G. Beneventano 

Departamento de Fi'sica and IFLP- Universidad Nacional de La Plata - 
CONICET 



Paola Giacconi 

Istituto Nazionale di Fisica Nucleare - Sezione di Bologna 

E. M. Santangelo 

Departamento de Fi'sica and IFLP- Universidad Nacional de La Plata - 
CONICET 

Roberto Soldati 

Dipartimento di Fisica - Universita di Bologna, Istituto Nazionale di Fisica 
Nucleare - Sezione di Bologna 



Abstract. 

We study 1-loop effects for massless Dirac fields in two spatial dimensions, 
coupled to homogeneous electromagnetic backgrounds, both at zero and at 
finite temperature and density. In the case of a purely magnetic field, we 
analyze the relationship between the invariance of the theory under large gauge 
transformations, the appearance of Chern-Simons terms and of different Berry's 
phases. In the case of a purely electric background field, we show that the effective 
Lagrangian is independent of the chemical potential and of the temperature. 
More interesting: we show that the minimal conductivity, as predicted by the 
quantum field theory, is the right multiple of the conductivity quantum and is, 
thus, consistent with the value measured for graphene, with no extra factor of it 
in the denominator. 



PACS numbers: ll.10.Wx, 02.30.Sa, 73.43.-f 



1. Introduction 

Planar electrodynamics, i.e. electrodynamics in two spatial dimensions, has attracted 
attention for years, mainly due to its theoretical departure from what we are used to 
in three spatial dimensions - for two reviews of the subject, with emphasis on Chern- 
Simons terms, and many relevant references see [TJ [5] . A pedagogical introduction to 
planar electrodynamics can be found in J3J. 

In turn, more than twenty years ago, it was shown that, in the tight-binding 
approach, the effective theory at low momenta for a two-dimensional lattice of carbon 
atoms (which is now known as graphene) is nothing but a theory of massless Dirac 
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fields in a .S+i-dimensional Minkowski space-time, in a reducible representation of the 
Clifford algebra (each irreducible representation corresponding to one Dirac point or 
valley) and with two fermion species, representing the spin of the interacting electrons 

H Ill- 
Stable monolayer samples of graphene were first isolated in 2004 [6J. The 
measurements, in 2005, of the Hall conductivity in those samples confirmed that 
massless planar quantum electrodynamics is adequate to describe the effective 
behavior of quasi particles in graphene [7]. The special, unique form for the 
conductivity of graphene, its many potential applications, as well as the link that 
graphene establishes between a relativistic field theory and a condensed matter system 
have triggered, during the last years, a huge amount of work, both experimental and 
theoretical, on the subject - for a recent review, see [8]. 

In [9] , we studied planar spinor electrodynamics at finite temperature and particle 
density, in the case of massless Dirac quantized fields in the presence of a homogeneous 
magnetic background. To this aim we evaluated, in the zeta function regularization 
approach [lOj , the one loop effective action of the three-dimensional Euclidean theory, 
at finite temperature and in the presence of a real chemical potential. From such 
a result, we obtained the mean particle density, or equivalently the net charge, as a 
function of the temperature and of the chemical potential. The zero temperature limit 
of this charge, followed by a Lorentz boost, allowed us to obtain the Hall current as a 
function of the chemical potential. In particular, the behavior of the Hall conductivity 
around zero chemical potential was shown, in the same reference, to depend on the 
selection made for the so called "phase" of the Dirac determinant [11] (for a discussion 
of this point, see also p~2] ). The present paper is a natural continuation of our 
previous work on the subject. Indeed, here we will study one-loop effective actions 
in the presence of those particular classical configurations characterizing the relevant 
experimental setups. 

In sections [2] and [3l we present some basic conventions and calculations needed to 
study the massless theory in a purely magnetic homogeneous background, this time for 
a general complex-valued chemical potential. The introduction of an imaginary part 
in the chemical potential will allow us to discuss, in sections 2] and O the invariance 
of the effective theory under "large" gauge transformations, their interpretation as 
4-7T rotations, and the relationship among different phases of the determinant, the 
appearance (or lack thereof) of geometric phases at finite temperature, as well as their 
interpretation as Chern-Simons terms or Berry's phases [13] in the zero temperature 
limit. These points were briefly discussed by two of us in the proceeding issue [T3] . 
Here, we present a more detailed analysis of these topics. 

Sections [6j [7] and [8] treat the complementary case of a purely electric constant 
background field, both at zero and non- vanishing, temperature and chemical potential. 
In section [6l we use the Bogoliubov transformation method to study the vacuum 
persistence probability, determined by the imaginary part of the effective action at zero 
temperature, in the 2+1 dimensional Minkowski space-time. We obtain the probability 
of vacuum decay and show that, as expected for massless particles, it is far from being 
small. In section we evaluate, through zeta function regularization, the effective 
action for the same background at finite temperature. We show that, in this case, the 
net charge is zero at all temperatures (equivalently, the effective action is independent 
of the chemical potential, whether real or complex). We also show that, in going back 
to Minkowski space-time, the imaginary part of the effective action coincides, at all 
temperatures, with the one obtained in the previous section in the zero temperature 
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case. 

The most impressive result in this paper is probably the one contained in 
section [5J where we evaluate, through the introduction of an adequate Lagrange 
multiplier, the mean current in the presence of a constant electric field. We show 
that such a calculation predicts a minimal conductivity which, at variance with many 
of the different results presented in the literature (see, for instance, [HI E]), is 
compatible with the measurements presented in [8J . Finally, section [9] contains our 
main conclusions, while the appendix presents a short overview of the main properties 
of the parabolic cylinder functions, needed to follow our calculation of the probability 
of vacuum decay at zero temperature. 

2. Dirac field interacting with a constant magnetic background at zero 
temperature 

In this section, we use the 2+1 Minkowski space-time metric = diag (+, — , — ) , 
and, throughout the paper (unless otherwise stated), we adopt natural units h = c = 
6=1, where c denotes the "speed of light" for the two dimensional system. Note that, 
for graphene, c is approximately two orders of magnitude smaller than the speed of 
light in vacuum [8j [17] . We choose the following representation for the Dirac matrices 

7° = (T3, 7 1 = io\ and — ia% . (1) 

They satisfy the required properties 

7 Mt = 7V7°, 7° = 7 0t , l k = ~l k \ k = 1,2; (2) 

{7 , ',7"}=7Y+7Y = VI- (3) 

As is well known [18] , a second, nonequivalent representation of the Clifford 
algebra can be chosen, for instance, by reversing the sign of one of the spatial Dirac 
matrices. We will always consider massless fields, two fermion species or flavors, and 
the direct sum of both nonequivalent representations, i.e., the effective model for 
graphene [H[S]. 

For the representation in equation ([1]), the 1-particle Dirac hamiltonian operator 
Hd is determined by the stationary solutions of the Dirac equation 

(^ + e4)*(x)=#>*(x) = 0, (4) 

where — e is the negative electron charge. From the previous equation, it is easily seen 
that 

H D = -eA a -a k D k , (5) 

where ^y°j k = a k — (— a 2 , oi) and and Dk = idk + eAk , k = 1, 2 . 

If a uniform magnetic field is present along the positive Oz— axis perpendicular 
to the Oxy— plane, after setting 

* £ (x°,x) = e- lEx ° i> E {x\x 2 ) = e~* m ^ E (x, y) , (6) 

and choosing the Landau gauge = (0, 0, — xB) , with B > , we get the 1-particlc 
Dirac hamiltonian 

Hd = id x a 2 - (id y + eBx) o\ 

(0 d x — idy — eBx 

—d x — idy — eBx 
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which allows us to determine the eigenvalues of the Hamiltonian from the set of 
equations 

(E — d x — p + eBx \ . , . 

d x - P + eB X E J =0. (8) 

In order to find the energy eigenstates, it is convenient to introduce the dimensionless 
coordinate 

£ = V2eB ( x - -?-) , d x = V2eBd i . (9) 
In terms of this variable, the 1-particlc Dirac hamiltonian can be written as 



H D = V2eB 
with 

s± = de t | 

where the raising and lowering operators satisfy 



[<5±, <5 T ] = ±1, [<5 T , <J T ] = 0, (10) 

8+\n) = -|n+l>, 5_|n)=n|n-l), (11) 

(z|n)=0„(z), (n'|n>=5 nn /. (12) 

Now, if we require the Dirac Hamiltonian to be self-adjoint on the domain L 2 (M. 2 ) , 
we are driven to the orthonormal set of functions 

(n\V2^y 1/2 D n (0 = MO, n = 0,1,2,... , (13) 

where 

D„(0 - 2-"/ 2 e-« 2 / 4 i/„(C/V2) (14) 

are parabolic cylinder functions of integer order. Notice that the normalization 
constants have been chosen to satisfy 



f 



d£ <l> m (0<Pn(0 = S mn , n,m = 0, 1,2, ... . (15) 



Hence, we obtain the following complete and orthonormal set of stationary 
improper eigenfunctions of the first order Dirac Hamiltonian. These comprise an 
infinite number of modes with zero energy, 

* o, P (t, x, y) = ip o, P (x) (Zero modes) 

1/4 
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and an infinite number of modes with nonvanishing energies, 
(t, x, y) = ^ (x) exp{± i E n t} 



eB\ 1/4 1 



exp{± i E n t + ipy} 
2tt 



(17) 



corresponding to eigenenergies ^fE n = ^f\ / 2neB , where tan 6^ = TE n . 

The spectrum is purely discrete and each eigenvalue exhibits the well known 
continuous Landau degeneracy per unit area, 

« 18 » 

In fact, the improper Dirac eigenspinors are normalized according to 

dx^J,,(x)V'o,fc(x) = 5(A;-0, (19) 

dx (x) 7 ° (x) = S mn 5(1 - k) , (20) 
and, as a consequence, we have 

dk V>o,fc(x)^o,fc(x) = A> 

dfc ^ ( x ) 7° Vfi (x) = A , V n E N . (21) 

The corresponding eigenfunctions in the other nonequivalent representation, 
where 7 1 — > — 7 1 , are obtained through -0 O , k ~^ ci^o.fci VVfc — * o'lV'n, ft, an d the 
corresponding energies are the same as the ones found in the present representation. 
Finally, to describe graphene, an overall degeneracy of 2 (two species or "flavors") 
must be considered to take the spin of the electrons in the original tight binding 
model into account. 



3. Dirac field interacting with a constant magnetic background, at finite 
temperature and density 

In order to consider the effect of temperature, we study the same problem as in the 
previous section, this time in the Euclidean three dimensional space, by considering the 
effective one-loop action in the path integral approach. This requires the evaluation 
of the determinant of the corresponding Euclidean Dirac operator, through some 
regularization method, which we will choose to be the zeta function one [10] . To 
this end, we take the Euclidean gamma matrices in one of the two nonequivalent 
representations as 

7j = -il' J , 3 =1,2, 7 3 = 7°, 

so that 

7>=0m> A* = 1,2, 3, {7^,7^1 = 2(5^1, 

and set 

ix = x 3 , x^ = (£1,2:2, £3) = (x, t) = (x,y,r), 
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A = iA a , = (A,A 3 ) = (0,Bx,-iv/e), 

where a generalized complex chemical potential v = [i + ij has been introduced. 
Note that the presence of a real chemical potential in the Minkowski space-time just 
amounts to a trivial overall shift of the energy scale. At finite temperature, instead, 
the derivative of the effective action with respect to the real chemical potential will give 
us the mean value of the net charge of the theory. The imaginary part of the chemical 
potential will, later on, allow us to interpret our results in terms of topological phases. 
To describe graphene we treat the massless Euclidean Dirac operator 

( P) Euci = °i ( _ * ®x) +< x 2(-idy- eBx ) + a 3 (-id T + iv) 

= <r\(p\ - eBx 5 2 \ - 7 8 3 \) + %ixu 3 . (22) 

In order to evaluate the functional determinant, we look for the spectral resolution 
of (|2"2"|) . which is not a self-adjoint operator for \x = 5ie v 7^ . We restrict its domain 
to regular square summable functions in x . Moreover, we shall adopt antiperiodic 
boundary conditions on the strip < r < f3 , where (3 = 1/T, (note that ks has 
been put to one in natural units). As a consequence of this antiperiodicity, the three 
dimensional Euclidean space manifold is the direct product M 3 = R 2 xC 1 (^r), where 
C n (r) is the circumference of radius r and integer winding number equal to n € Z. 
Hence, the symmetry group of .M3 becomes the direct product U (1) x 0(2) = U(1) x 
U(l) , that is the square of the unitary, non simply connected, Abelian group U(l). 

To satisfy the antiperiodicity, we introduce the Matsubara angular frequencies 

ui = (2t+l)jj, t€Z, 



and propose, for the eigenfunctions of (|22 

1 

7w 



^(x,V,t) = Y] exp{i(TU)i+py)} %j) tp {x), (23) 



where 

dp . 



E 

— — 00 



) ( z ) 



= 



Thus, the problem to be solved becomes 

(u>i + iv — A — i d x + i ( eBx — p] 

— id x — i(eBx — p) —oj£ — iv — \ 

The resulting spectrum is complex and consists of two parts. We will call the first part 
asymmetric, since, given one eigenvalue A in this part of the spectrum, —A doesn't 
belong to the spectrum. The second part, which we will call the symmetric part of 
the spectrum, instead, is such that, for each eigenvalue A belonging to this part of 
the spectrum, —A also belongs to the spectrum. ( Note that this last part behaves as 
expected from the eigenvalues of a square root operator). 

(i) The asymmetric piece of the spectrum, which is an infinite "tower" of states 
associated to the lowest Landau level in the Hamiltonian, is given by 

A^ = cue + iv = uit — 7 +ifi, (leZ) (24) 
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and has a corresponding set of eigenfunctions 

T , , feBY ( Q-keB(x-p/eB) 2 \ 

*, P ( W )=(-j [ j 

_ 1 

x cx.p{iru;i + ipy} (27r/3) 2 . (25) 
Their normalization is given by 

r- ft rOO pOC 

dr dx dy^ ei p(x,y,T)^* q (x,y,T) = 5 ej S(p-q), 

JO J — oo J — oo 

and the degeneracy per unit area is, as in the zero temperature case, 

■ ft f-OO 



J drj dp \^ iiP (x,y,T) | 2 = — = A . 



(ii) The symmetric part of the spectrum is related to the excited Landau levels. It is 
given by 



A = ±A/,„, A/, n = y \\ + 2eBn, (26) 

with n € N and I e Z . To obtain the corresponding eigenstates, it is convenient 
to set 



z = xV2eB - p\ — , d x = V2eBd z , (27) 
V eB 

so that 

d x ± (p- eBx) = V2lB (d z T |) = V2^B 5± , (28) 

where the raising-lowering operators satisfy 

[S±,S T ] = ±1, [S T ,S T ] = 0, (29) 

6+\n) =-\n+l), S-\n) =n\n-l), (30) 

(z\n) = <p n (z) , (m\n) =6 mn , (31) 

_ 1_ 

<j> n {z) = (n!v^7r) 2 £>n(2), n = l,2,.... 
Hence, if we rewrite the Euclidean Dirac operator in the form 

(W Euc i = ^l f ^ 1 , (32) 



where = Xg/w2eB , it turns out that the normalized eigenstates are given 
by 

Vt, P ,n,±= (2eB)i exp{iTUe+ipy} (4tt/3)~5 

x (l + |x,,„, ± | 2 )-^ f ±i4> 2 {Z) ( , } . (33) 

with 
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( A ^ „ T A < 



n V2eB 
A| jri ± A^ 



Moreover, by setting 



one finds 



and, thereby 



A^,„ = p^,n (cosa^„ +i sina^,„) 



(c^-7) -^ 2 + 2eBn + A^i 2 (a; e - j)' 
2^(0)^-7) 



tan 2a 



3?e X£ 



n,± 



(w^— 7) —/i 2 + 2eBn 



pl, n cosat, n =F (w* - 7) 



^mx^ ni± = /) (j „sina (jn T A 4 ■ 

It might be convenient to introduce another different parametrization for the 
eigenspinors, which makes use of a pair of sequences of angles, namely, 



( i e z , n e N ) 



Thus, after calling 



we can also write 



xi,n,± = tan0^ n exp{i9^ n }, 



^£, P ,n= ^ 47r/j exp{zrw £ + zpy} 
±i 0„(z) cos9f n 



exp{iO e n } 4> n -i(z) 



smt 



(34) 



Also in this case, the degeneracy per unit area is given by the Landau factor A 
in psp . We remark that, when going to the other non-equivalent representation, the 
symmetric part of the spectrum remains the same, but the asymmetric part transforms 
according to A^ — » —X( or, equivalently, v — > —v. 

4. The effective action. Invariance under large gauge transformations in 
the case of a complex chemical potential 



It is important to note that, as can be seen from (|24[) and (|26[) , the whole spectrum is, 
in any of the two representations of the Clifford algebra, invariant under the so called 
large gauge transformations [TJ 119) 

v -> u' = u + ^H V(x) -> V'W = V'(x) e lr ^ , (j£Z) 
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where x = (x, y, r) , which are precisely those Euclidean gauge transformations leaving 
untouched the antiperiodic boundary conditions satisfied by the spinor Euclidean 
wave functions and, in consequence, the Fermi-Dirac statistics. In other words, the 
discrete imaginary shift of the chemical potential j3v — > (3v + 2irij (j £ Z) is 
a symmetry transformation of the theory. On the contrary, any other imaginary 
shift of the chemical potential cannot be compensated by a phase tranformation of 
the Euclidean spinor wavefunction, for the latter one would spoil the antiperiodic 
boundary conditions that the Euclidean spinor field must satisfy in order to have the 
right statistics. 

The Euclidean effective action for a sample of area £1 , as defined in terms of the 
zeta function regularization technique [lOj , is given by 



T cS = lnZ = \ndet(l (]p) 



Eucl 



as 



s = 



= - C(o,pe ) 

= ClA[A(a) + S(a,b)] . (35) 
In the last equation, we used the dimensionless reduced variables 

the arbitrary length scale £q being introduced to render the argument of the 
zeta function dimensionless. It is important to stress that physically meaningful 
results must be independent from £q. Here A and S denote the dimensionless 
reduced contributions to the full effective action that originate respectively from the 
asymmetric and from the symmetric part of the spectrum of the Euclidean Dirac 
operator. 

In reference [9] , we presented a detailed calculation of the effective action in the 
case of a purely real chemical potential (v = fi), and showed that different selections 
of the cut in the complex plane of the eigenvalues of (|22|) during the evaluation of the 
determinant lead to predictions for the Hall conductivity which reproduce the behavior 
measured for mono- and bilayer samples. In this section, we will concentrate on the 
contribution to the effective action due to the asymmetric portion of the spectrum in 
the case of a complex chemical potential- a survey of this calculation appears in |14j . 

The contribution to the effective action coming from the symmetric part of the 
spectrum does not suffer from regularization ambiguities. It is given by 

oo 

S(a, b) = 2b (r H) + E ln {(! + ze ~ 0Sn ) i 1 + z- 1 ^)} , 

n=l 



where z = e 2va , e n = \/2eB n . and £r is Riemann's zeta function. Notice that this 
expression is invariant under the discrete symmetry v — > — v, which is a symmetry of 
this portion of the spectrum. 

On the other hand, the contribution to the effective action arising from the 
asymmetric part of the spectrum is given by 



A(o) = -Ca(Q,o) . 
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with 



Ca(s, a) 



( 2n£ 
\ 



-s oo 
-s oo 

£ 

1 = 



Bv 



1 Bv 
- + i — 

2 2tt 



{i + q+y s + cT™ s (t + q _y 



(36) 



The right hand side involves the sum of two series representations of Hurwitz zeta 
functions, £(s, q) , which are well defined iff 



In 



5ies > 1 , Bj ^ ±?r(2n- 1) , 

we have introduced the complex quantities 

iB 



1 

Z Z7T 



IT < Bj < IT 



The regularization ambiguity, which gives rise to a phase of the determinant |llj , 
has its roots in the factor e Tl7TS or equivalcntly, in the selection of the cut of the 
complex power in the plane of the eigenvalues of (|2"2"|) . Selecting one or the other sign 
(cut in the lower(upper) complex plane), irrespective of the sign of /x, we get 



(a(s, a) 



2ir£ { 



t =p 17TS 



C(*,«-)} 



For the contribution to the effective action we obtain, after evaluating the derivative, 
A (a) = ln27r- ]nT(q + )T(q_) ±iir (I - q_) 

= In 2?r - In T (| + ia) F (| - ia) T vra , (37) 
for -7r</37<7r, fieM. 

The previous result can be also be cast in the form 
A{a) = In 2 cosh na =F ira 



In 



e 2 



/3( M +i7) 



T ^ Z 3 ( A» + H ) 



(38) 



So, in our case, a well known sign ambiguity [111119] appears when evaluating the 
determinant of the Dirac operator, as it happens whenever this last operator is not 
the square root of a second order operator of the Laplace type. Note that, if the so 
called "phase of the determinant" were chosen with the same sign for all values of /j,, 
then the contribution of each representation to the effective action would not be an 
even function of fi for 7 = 0. As a consequence, the contribution to the charge due to 
each Dirac point would not change sign as [i — > —fj, . 

Now, it is always possible and convenient (although not mandatory) to relate 
the two opposite phase prescriptions with the sign of the (real part of) the chemical 
potential. 

Turning back to eq. (|36|) , if we cut the complex plane of the eigenvalues A f = 
(£ + q±) (£ = 0,1,2,...) along the upper half plane for positive values of \i and 
along the lower half plane for negative values of /i, for example, we obtain what we 
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will call the standard prescription, and will identify through the index k 
this case we have 



7r a . 
7r a . 



A(a) — In 2 cosh ira = 
We shall call the non-standard prescription, k = 
A(a) — In 2 cosh na = 



+ ir a . 

— 7r a . 



for fi > , 
for fi < ^ K 

-fl, the opposite one 

for fi > 
for n < 



1 



1. In 



(39) 



( k = +1 ) . (40) 



In so doing, for 7 = 0, we recover our previous result [9], viz., 



r eff = A {| k/3| M | + In 2cosh ^ + 26 C (-5) } 

OO 

+ A ^ In { (1 + z e^ £ ") (1 + z- x e-^) } . 



(41) 



These two criteria lead, for a complex chemical potential such that —tt < /?7 < tt , 
to a contribution of the asymmetric part of the spectrum to the effective action, which 
is given by 

A(a)= In ( e ™ + e" ™ ) + nan sgn(/z) 



= In 



+ |/ t /3(/i + i7)sgn(/i) . (42) 



This contribution turns out to be invariant under the generalized conjugation 
symmetry a — > — a. 

One can take the limits /?7 — ► ± 7r in the last expression and write 



lim A(a) = In 2sinh 

/?7 — > 7T 2 



lim A(a) 

f3~f — > — 7T + 

and therefore 



In 2 sinh 



1 „ , vr't , 

Ik/3\iil\ + -(1 



y(l 



K ) sgn(ju) , 
h k ) sgn(//) , 



lim — lim 



A(a) = 



for k = -1 

2vTi sgn( /1) for k = +1 



However, the exact calculation can be performed at both extreme points of the 
interval, and both results coincide, as dictated by the invariance of the spectrum under 
large gauge transformations, which amounts to f3"f — > (3j + 2n. The zeta function 
regularization, which is based solely on the properties of the spectrum, respects this 
invariance for all values of pj. In a zeta function regularization approach, not only 
the partition function is invariant under these transformations: it is the effective 
action itself that remains invariant. So, we can limit ourselves to give the result for 
the contribution to the effective action of the asymmetric part of the spectrum, for 

— 7T < f3j < 7T , 



Q^Teff = A [A (a) + S(a,b)] 
pR 

= — {ln(e™ + e~™) + 7ra/csgn(/i) + 2&£ (-f) 

OO ^ 

+ ^ln[(l + ze-' 3£ ")(l + z- 1 e-' 3£ ")] . (43) 
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The result for any other value can be obtained by using 

T eS (Pn + ipY)=T eS (Pn + 2iri3 + il3i) , (j6Z). 

In fact, for any of the two selections of k, the whole effective action, per unit degree 
of freedom is invariant under large gauge transformations. As explained in [14] , this 
amounts to an invariance under rotations of angle 4-7T around an axis perpendicular to 
the plane of the sample in Euclidean space. We discuss this point in the next section. 



5. Zero temperature limit 

Owing to the fact that the eigenfunctions corresponding to the asymmetric part of the 
spectrum are eigenstates of 03 with eigenvalue +1 , one can equivalently write gauge 
transformations in the form 

* £, P (x,y,T) -> exp {^ia 3 (27) r} (x,y,r) . 

This last expression shows that, as the Euclidean time r grows from to /3, spinors 
are rotated by 2/37, since i 03 is the generator of rotations in the Oxy~ plane. In 
particular, 0^ — 2tt corresponds to a rotation of angle 4-7T around the magnetic field. 
Such a rotation must be an invariance for spinors, which provides one interpretation 
of the invariance under large gauge transformations. 

On the other hand, (3j — tt corresponds to a rotation of angle 2ir . As can be 
seen from (|43[) . at finite temperature, such a transformation turns the statistics from 
Fermi-Dirac to Bose- Einstein, for any of the two selections of phases in one of the two 
inequivalent representations (for a discussion about Bose-Fermi transmutation and 
relativistic "anyons" see, for instance, [5U]. For k = +1, it also gives rise to an overall 
phase of tt per unit degeneracy in the partition function, which can be recognized as a 
Berry's phase. Such a phase is the contribution which survives in the low temperature 
limit. In fact, in this limit one has from equation (|4"2"|) , 



nAA{a) = i(K + l)ftA/3(^ + i7)sgn(», (P~>oo) .(44) 

Thus, in the zero temperature limit, the effective action vanishes for the standard 
phase prescription (k = —1), so that the partition function is rotationally invariant. 
For k = +1, the imaginary term in QAA(a) (i.e., the one proportional to 7) is the 
Abelian Chern-Simons term 

if3 1 n— sgn( M ) = .4 CS = k cs — e^F^A,, , 

Z7T Z Z7T 

where kcs denotes the Chern-Simons coefficient, a pure number of a topological origin 
(remember that ie^ = fj, + vy). Our calculation clearly shows that the latter one is 
invariant, as expected, under large gauge transformations, which represent rotations 
by integer multiples of An . On the contrary, under rotations of angle 2ir (f3j — tt) the 
partition function is multiplied by (when taking into account the two flavors) 

Z(/? 7 = tt) Z e ±2 ™ OA . 

Hence, the requirement of the partition function, in each representation, being 
invariant under rotations of angle 27r, for an area f2, entails a quantization condition 
on the coefficient of the Chern-Simons term or, equivalently, on the reduced flux of the 
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magnetic field, i.e., HA = N e N. Interestingly enough, this is precisely the condition 
for physical states to transform according to the one dimensional ray representations 
of the magnetic translation group [3T]. 

Note also that the geometric phase appearing for k = +1 leads, in the zero 
temperature limit, to a Berry's phase which equals 7r for each degree of freedom. 
A similar calculation in the other irreducible representation of the Clifford algebra 
would give exactly the same results, if the same criterion were used to fix the phase 
of the determinant. As stressed in [9], the right behavior for the Hall conductivity in 
monolayer graphene appears when choosing opposite criteria in the two noncquivalcnt 
representations, which corresponds to a total Berry's phase of 7r per unit degree 
of freedom. On the other hand, using the same criterium (k = +1) in both 
representations, which corresponds to a Berry's phase of 2ir per unit degree of 
freedom, correctly reproduces the Hall conductivity of bilayer graphene [9j. For 
other discussions of Berry'a phases in mono- and bilayer graphene see, for instance, 

El M M- 



6. Spinors in a constant electric field : Pair production at T = 

Consider now a graphene sample in the presence of a constant electrostatic field 
pointing towards the positive Ox— axis, i.e. F 10 = Foi = E x = E > . In this section, 
we will choose a gauge which leads to non-stationary sets of solutions. Throughout 
this section we will present a calculation, for the 2+i-dimensional case, which is quite 
similar to the one performed in [25] for the 3-/-i-dimensional Minkowski space-time. In 
particular, to solve the Dirac equation in the present ^-/--/-dimensional massless case 
it is convenient to employ the following representation for the gamma matrices, 

7° = cti , 7 1 = ia 2 , 7 2 = i<J 3 ■ (45) 
After setting (x°, a; 1 , x 2 ) = (t,x,y) we get the massless Dirac operator in the gauge 
[25] A» = (0,-M,0), i.e., 

Tj) = id t o\ + (id x + eEt) ia 2 - d y cr 3 , (46) 

which is explicitly time dependent. In order to obtain the solutions in this gauge, it 
is convenient to introduce the partial Fourier transforms 

-i />oo P OO 

*(*, x,y) = — dp dk e ipx+2ky <f(t,p, k) , (47) 



2?r _ 

with p = (p, k) as well as the dimensionless quantities 

^=(p-eEt)/y/eE, X = k 2 /eE. (48) 

If we write 

•<«•*»>- M)- »ft»*>- ( !;^) ) • « 

(50) 



wc obtain the coupled differential equations 

ik <p + VeE (id$ + £) x = , 
VeE (-id£ + £) <p + ik x = . 
Then, we can write 

X=(*/>/A) + ( A ^°) ( 51 ) 

(d| + £ 2 + A + t) </? = . (52) 

(ide + 0x = o, Hde + 0£ = o, (A = o). (53) 
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The general solutions of eq. (|53p can be expressed in terms of two different sets 
of linearly independent solutions, involving parabolic cylinder functions |26j . 

A) = A ia D_ iX/2 [+ (1 + i) f ] + D tt/J _, [+ (1 - i) £] , (54) 

where Ai n and i?i n are complex constants. From the recursive relations 

(-id £ + £>_ iA/2 [ ± (l + i) £ ] = T ^ (l + *) D_ iX /2-i [ ± (l + ») £ ] , 

+ -Da/2-i [ ± (1 - i) e] = ± (1 + *) £>iA/2 [ ± (1 - £ ] , (55) 

we get the normalized spinors, which are solutions of the original massless Dirac 
equation and read 

1 f . ttA ' 



(56) 



-* p (i,x) = — exp -j ip • x - — 



2Z?_ iA/a [(l + i)£] 
+ (l-<)>/Ai3-a/2-i[(l + *)^] 



(57) 



. , i r. tta 

+ *p (*,x) = — exp | ip ■ x - — 

-(l + *)VAI>iA/2-l[(l-<)e] 

2Da/2[(l-i)f] 

An alternative set of linearly independent solutions can be selected, that is 

A) = A out D^ x/2 [- (1 + i) £] + 5 D * ut Dix/2-i [ - (1 - *) £] . (58) 
Then, the recursive relations in equation (|55[) give 

i T . . 1 f 7rA 

*p (*, x) = — exp j«p • x - — 



( 



2£>_ iA /2 [-(! + »)£] 

(i- J )VA J D_ 4A/2 _ 1 [-(i + l )e] 

1 f ttA^ 
*p (*' x ) = "TZ ex P i 1 P ' x ~ ~ 



(59) 



4tt 



(60) 



r +(i+i)V\D iX/2 ^[-(i-i)t} ) 

{ 2D iX/2 [-{l-i)i] J • 

The above spinor solutions fulfill the orthonormality relations 

(±* P , ±* q ) = / ±*p (*« x ) ±*q (*> x ) dx = 5 (p - <*) > ( 61 ) 

p , ±* q ) = y ±* t ( t) x ) ± Wq ( t) x ) dx = 5 ( p _ q ) , ( 62 ) 

( ± *p, T * q )=0=( ± * p ,^ q ) . (63) 
Note that for A = = k , owing to Dq(z) — exp{ — z 2 /4} , both pairs of solutions 



* p (t,x)= +*p(*,x) = i- cxpizpx- U 2 (t)\ fil, (64) 



coincide and read 

* p (t,x) = -* P (t,x) = cxp|zpx+ ^ 2 W| ( J ] ■ (65) 
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By looking at the asymptotic behavior of the parabolic cylinder functions, one finds 



-*p (*,x) 
+ * P (t,x) 



1 

2^ 
1 

2^ 



2 f2(t)]-*V4 

2£ 2 (i)] jA/4 



+* p (t,x)~-[2e 2 (t)]- a / 4 

Z7T 



(i) 
(?) 

m 



1«P-x--C (*) 



cxp<{ tp x- 



(*,x) 



^[2£ 2 (i)p A / 4 J I exp^p-x+-f(/i 



(i->-oo) (66) 

(t-v-oo) (67) 

(t->+oo) (68) 

(i^+oo) (69) 



Hence, we understand + v f , P as the wave function of an incoming quasi-particle, while 
we associate the wave function to an incoming anti-quasi-particle. Conversely, 
we shall associate and "^p to the corresponding outgoing wave functions [2"5~] . 

The above sets of incoming and outgoing solutions are related throughout a 
Bogoliubov-like unitary transformation, that is 



We have 



where 



with 



+ ^ p {x) = ci + ^ P {x) + c 2 ~^ P {x) 
■*pW = -c* 2 + V P (x) + cj -*p(a?) 



C2 



(+* p , +* q ) =JVj<5(p-q), 

iVA^(p-q), 
exp{— ttA/2} 5(p — q) 



"tfp, _* q ) = JV A <J(p-q), 



/A ^ fi\\ . , /vrA\ f ttA 3™ 



N x \ 2 



1 — exp {— 7rA} 



(70) 

(71) 
(72) 
(73) 

(74) 



Then, using the above listed orthonormality relations (|6ip . ([62]) and (|63|) . we 
immediately find 

ci = iV A * , c 2 - exp {- ttA/2} = c* . (75) 

As a consequence, in the transition to the quantum field theory, we can write the most 
general operator solutions of the Dirac equation as 



tfinfs) = dp [ai„(p) +^ P (x) + b\ n (p) -* p {x) 



*out(«) 



dp 



flout (P) + *p(a:) + Oout(P) *p(«) 



where each set of creation and destruction operators fulfills the standard canonical 
anticommutation relations. From the condition ^^(x) = ^ ou t{x) one has 
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flout (p) =c 1 a- m (p) -c 2 6 1 t 1 (p), (76) 

6L(p) = ^om(p) + ^6to(P). ( 77 ) 

so that 

<0pout| = (0pm|6 in (p)6Lt(p') 
= (0 pm|a in (p)ai ut (p') 

= (0pin|(c^(p-p')-C2ain(p')6in(p)) • (78) 

Once we have the properly defined asymptotic in and out spinor states, which are 
explicitly gauge dependent, we can naturally define and calculate the gauge invariant 
probabilities for the processes of pair creation and annihilation in the presence of 
the background electrostatic field. We start by obtaining the vacuum persistence 
amplitude at a given p. From equations (|75|) and (|75|) we find 



(0pout|0 qin) = c[ S(p — q) 

= Jdx-%(x)-*^x) = N x 6(p-<i), (79) 

with N\ given by equation ([74]) . 

In turn, the probability of vacuum decay at a given p is 

w x = 1- \N X \ 2 , 

or, more explicitly, 

k 2 

w x = exp{- ttA} , A = — — . (80) 

eh 

It is clear that the latter one is also the probability of the inverse process, i.e. the 
annihilation in the state p with transfer of energy to the external field. As a matter 
of fact, we also have 

j dx+tft (t,x)_* q (t,x) = *(p-q) cxp{-i^A} . (81) 

The total vacuum persistence amplitude can be written, after recovering physical units, 
in the form - see the recent up to date review [27], 

(out0|0in) = expj^ [%teT cS {E) + i3mr cff (£) ] j , (82) 

and the corresponding probability is 

| ( out | in ) | 2 =expj — QmC eS (E) \ , (83) 

where C c s(E) is the effective Lagrange density in the presence of a background 
electrostatic field E , whereas T is the total time, and fl denotes the area of the sample. 
It turns out that T e g(E) contains a real part that describes dispersive effects like the 
vacuum birefringence, as well as an imaginary part that concerns absortive effects 
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like vacuum decay. Now, according to the above derivation, the vacuum persistence 
probability, can be obtained from 

(fiT)- 1 ln|(outO|Oin)|^ i ^|^fcln(l-exp|-^|}) (84) 

where we have introduced the number of states between p and p + dp, which is , 5L 2 , 
and we have used that dp = e E dt. 

So, for the imaginary part of the effective action, we have, after taking the 
degeneracy / into account, 

l^L eS (E)=-f^P.V.J_ 




where c is the speed of light for the two dimensional graphene sample and / is the 
number of fermion species or flavours, which is four for graphene, since the result is 
the same in both nonequivalent representations of the Clifford algebra. Moreover, the 
symbol P.V. stands for the Cauchy principal value prescription, i.e., P.V. J_ dk = 
2 lim e ^o / e °° dk , which is the natural even prescription for the integrand around k — 0. 
For previous, related results, see [28],l2"9]. 

Note that the vacuum persistence probability exp{(— 2/h) ft T3m£ c (f(£)} is 
quite small, even at finite time and for a finite area, due to the massless character 
of the quasi-particles in graphene. So, the probability of the complementary event 



P= 1- cxp|^- -nT%m£ cS (E) 

is high and by no means approximatively given by the effective action as in the case 
of massive particles/antiparticles pairs. The vacuum persistence probability obtained 
by replacing (|85[) into this last equation coincides with the zero-mass limit of the one 
reported in [30] ■ 

Some authors [3T], starting with J. Schwinger [35], call this quantity the 
probability of pair creation. Others (see [30J [33J and references therein), following 
A.I. Nikishov 34J, prefer to characterize the creation of pairs through the rate of pair 
production. Note that the rate of pair production was evaluated in the same number 
of dimensions, for the massive case, in [30]. In the zero mass limit, the result found 
in this last reference is given by the first term of the series in |85|) . This is exactly 
the result one would get after integrating (|80|) over all impulses, when counting the 
number of states between p and p + dp as before. 

Apart from the way one characterizes vacuum decay, the important point is, 
indeed, that the probability of actually detecting vacuum decay effects is another 
appealing characteristic of graphene, at variance with the case of massive particles, in 
which vacuum decay is very hard to be unraveled (see |27]V 



7. Dirac field interacting with a constant electric background at finite 
temperature and density 

Consider now the problem of a graphene sample in the presence of a constant electric 
field pointing towards the positive Occ-axis, i.e. F 10 = Fqi = E x = E > , at finite 
temperature and density. 
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As done in section^ here we will employ path integral methods and zeta function 
regularization. In this section, at variance with what has been done in the previous 
one, we choose a gauge which leads to a stationary set of solutions. 

To solve for the eigenvalues of the Dirac equation in three dimensional Euclidean 
space, with metric (+,+,+), it is convenient to employ the following representation 
for the Euclidean gamma matrices, 

7i = °i i 7 2 = CT 3 , 73 = G i ■ ( 86 ) 

After setting ixo =t, x^ = (xi, X2, X3) = (x,y, r) and taking into account that 
A^ = (0, 0, £x + iv/e) , (e > 0) with A Q — iA% , E = i £ , we get the Euclidean Dirac 
operator 

-^Euci = (id T - e£x - iv) er 2 + id x o\ + id y 0-3 . (87) 

In order to evaluate the effective action at finite temperature, one must impose 
antiperiodic boundary conditions on the eigenfunctions, as in section |3l In this gauge, 
it is convenient to write the partial Fourier transform 



1 r°° 
9(x, y, t) = —— V exp{i™,} / dk e ik * Hf t (x, k) , (88) 

where, as before, 




denote the Matsubara angular frequencies. 

Moreover, it turns out to be very convenient to introduce the dimensionless 
quantities 



& = (X e + e£x)/V^£ , d e =4r (teZ), (89) 

with \i — toi + iv , so that the eigenvalue equation reads (note that, here, A represents 
the eigenvalues of the Euclidean Dirac operator, and is different from A in the previous 
section) 

^ f W/T I *<(*,*) = A *<(*,*). (90) 

^ id,i -lit k/Ve£ ) 

If we write once again 

■ ( g&| ) ■ w 

we are lead to the coupled differential equations 



e£ (d e + xi + i (A + k) tp e = , 
e£ (dt — (H)tpt + i(\— k)xt = . 



(92) 



After solving these equations, 
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(i) When A = k , we obtain the following infinite set of degenerate normalized 
improper eigenspinors : namely, 



*„.,„, {.,:„. r) = ( — ) <'x.|){ /YvC, - //.•// )• I 2;< /) 



exp | — | e£ (x + Xi/e£ ) 2 | 



which fulfill 

/>oo />oo 

dr / dx / dy *J (x,y,r)^o t k,i(x,y,T) = 8i m 8(k-p) 



oo J — oo 



(93) 
(94) 



Their degeneracy is 

oo 

X! *o,fc,f ( x ^' r ) *o,fc,£ (^,?/,t) 



1 

2tt/3 



— 7 . ex P 



! = -oo 




where we have suitably introduced the rescaled and dimensionless complex coordinate 
z = (e£x + iv)J^ , while use has been made of the Euler-McLaurin summation formula 
- see for example equation (3.6.28) in ref. [3^1 ■ Note that I— th term of the series is 
an even function of I. Thus, all the odd derivatives vanish at I = 0. Moreover, the 
the (.— th term and all its derivatives vanish at infinity. As a consequence we come to 
the result that the number of degenerate eigenspinors per unit Euclidean "time" and 
unit x- length is provided by As = e£ /An 2 , the degeneracy factor being the same for 
all the remaining modes. 



(ii) For \ k we can write 
„ iye£ 



X 



(dj 



\—k 
-& 2 + A-l)^ = 0, 



A 



A 2 - fc 2 
e£ 



(96) 
(97) 



As in section [31 the eigenvalue problem for each £ £ Z is the one of a linear harmonic 
oscillator, so that the spectrum takes the symmetric form 



A fe „ = ± \Jk 2 + 2e£r 



n = 1, . 



.,00. 



k e 



(ye e z). 



The somewhat surprising conclusion we can draw from the above analysis is that 
neither the eigenvalues nor the degeneracy depend at all upon the temperature and/or 
the chemical potential. Hence, once we have at hand the spectrum of the Euclidean 
Dirac operator and its overall degeneracy, which is the same for all eigenvalues, we 
can turn to the evaluation of the Euclidean effective action, within the zeta function 
approach. In this case we have 



C( s ,ip£ y ) = A £ pe x e y 



(e y k)- s + 



i + C T" s 



\ ^ 



(2e£l y 2 ) 



,/2 

71=1 



2e£ J 



(98) 
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Here, the length scale t y has been introduced to render the second argument in 
the zeta function dimensionless. 

Now, eq. (|55|) can be rewritten as 



dk (£ v k)- s + (99) 



r(s/2) (2e££ y 2 ) s/2 



/oo poo - v 

dk / dt t*' 2 - 1 c ~ 
Jq n=l 



tn-tk 2 /2e£ 



For dies > 1 the two integrals in the second line can be interchanged and, if the integral 
in the first line is understood in terms of the Cauchy principal value prescription, 
according to our comment in section 21 we get 



- 1 \ / s - 1 \ 1 

(100) 



e^O 1-s T(l + s/2) v ' V 2 / 

We finally obtain the Euclidean effective action per unit volume in the form 
(note, however, that the contribution coming from the first term was defined here 
as — lim e ^o / c °° dk-^\ s _ i^y i the integrand being evaluated at 3?es > 1 and 
then extended to s = before taking the principal value limit) 



1 



ds 



s = 



where (see [26] ) 



ir^>" "~\ A J nJ2 

n—1 



1 oo 



As a consequence, we can recast the Euclidean effective Lagrangean as 

C " = W2~ (6f)3/2 ' (101) 
so that, after turning back to the 2+1 dimensional Minkowski space-time by means 
of the replacement £ = — iE , we obtain the effective Lagrange density 

C™(E) = <M(eEf/ 2 {(l + i)}. (102) 
It follows that, as in eq. (|83p . we end up with the correspondence 



which leads, after recovering physical units, to eq. (|85|) . i.e., the 2+1 dimensional 
counterpart of the celebrated Schwinger formula [53] 

a "^»< E > = '^Sw- (104) 
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where c is the speed of light for the two dimensional graphene sample and / is the 
number of fermion species or flavours, which is four for graphene, in perfect agreement 
with our previous result (|85[) . Moreover we find that in the 2+1 dimensional case the 
following remarkable equalities hold true, viz., 

Be/$(25) = 3m£ e M ff (£) = - £ c E ff (£) . 

8. Minimal quantum conductivity of graphene from planar QED 

Of particular interest in the presence of a constant electric background is the mean 
current density or, equivalently, the minimal conductivity of graphene. In order to 
evaluate the mean value of such quantum current density, we go back to our calculation 
in the previous section, this time introducing a Lagrange multiplier a through the 
non-trivial replacement k — > (fc + ia) . The derivative of the effective action with 
respect to the external parameter a will then give us the required mean value of the 
quantum current density. Actually, we can say that the control parameter a plays a 
very similar role to the one played by the chemical potential in the evaluation of the 
particle number or charge density. 

If we consider the manifold to be compact in the y direction, all we have to do is to 
turn back to equation (|100[) . where we have, after a shift in the transverse momentum 
integration variable k , 

C{s,p£ y ;a) = A £ ni3£ y ~ s \ [ dk(k + ia)- s + (105) 



l + cT™ 



r(s/2) (2e£) 



</2 



/OO p OC - x - 

dk / dt t 8 ' 2 - 1 V c 
-oo Jo n=1 



- tn-t(k+ia) 2 /2eS 



Now, the shift in the second term is irrelevant, since the Gaussian function is 
an entire function. Then, the a— dependence is solely due to the first term and, 
consequently, the only contribution to the minimal quantum current density will arise 
from the quantity 

/oo 
dk (k + ia)- s 
-OO 
P OO 

= A £ n (3i y - s / dk [(k + ia)- s + {-l)- s (k~iay s ] . 
Jo 

At variance with the case of the previous section (a — 0), the zeta regularization 
of this term is now well defined, the only ambiguity being the selection of the phase, 
much as in the calculation of the charge when the sample is subjected to a magnetic 
field. As in that case, one can make two selections of the phase. Indeed, 

poo 

Z s {s;a)=A £ n/3i y - s dk (k + ia)~ s + e -™ Ksi g"(«)(fc _ i a )- s ,(106) 



jo J 

where, as before, k = — 1 represents the "usual" selection of phase, whilst k = 1 the 
opposite "unusual" one. Thus, for 3?es > 1 , we get 

Z £ (s ; a) = -A £ Q [l - e " 2 ™ SK sign ^l . (107) 

1 — s L J 
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Performing the s— derivative and evaluating at s = 0, we get for the a— dependent 
part of the Euclidean effective action, 



rjf(a) = -2TTA £ np K \a\ = -^-/30«;|a| . (108) 

To the aim of recovering the minimal quantum current density it is necessary, first, to 
go back to the 2+1 dimensional Minkowski space-time, i.e., 

p F 

T™(a)= — Tnn\a\. (109) 

Now, it turns out that performing the derivative with respect to a , dividing by f2 T 
and multiplying by the elementary charge, — e, we have 

e 2 E 

(J m in)=--^—Ksign(a), (110) 

Note that, in the other representation, the result is the same if the same criterium 
is chosen to define the phase of the determinant. Thus, summing up the contributions 
from both representations, contrary to the case of the Hall conductivity, with the same 
criterium, multiplying by the two spins (flavors), and recovering the physical units, 
we find for graphene 

4e 2 

( Jmin ) = - -r- E k sign(a) . (Ill) 
h 

The quantization of the minimal quantum conductivity 

4e 2 

own = j- Ksign(a) (112) 

in terms of the quantum unit of conductivity makes this prediction entirely different 
from the results obtained, for instance, through the Kubo formula [T5l |T6] , where an 
extra factor of ir in the denominator appears. Moreover, for k = — 1 and positive a one 
obtains exactly the result in [S]. It's interesting to note that a similar experimental 
value is found also for bilayer graphene [|24 . 

As already stressed, this result is independent both from the temperature and 
the chemical potential. 



9. Conclusions 



In this paper we have carefully analyzed the response of a graphene sample under the 
influence of homogeneous electric and magnetic fields. Our present investigation has 
its roots in the well established quantum field theoretic model of a free massless Dirac 
spinor in two space and one time dimensions, i.e. the massless planar spinor quantum 
electrodynamics. In spite of its simplicity, this quantum field theory does exhibit 
remarkable features. As a matter of fact, in the presence of a uniform magnetic 
field at finite temperature and density, a nice connection has been established, in 
this manuscript, among the phase prescription of the effective action, the large 
gauge transformation invariance, Berry's phases and Chern-Simons topological terms, 
besides the different forms for the Hall conductivity we had already elucidated in our 
previous paper [9]. 

Moreover, in the presence of a constant electrostatic field, a derivation of the 
imaginary part of the effective Lagrangian at zero temperature has been presented. 
Our result coincides with the zero mass limit of the one obtained, for instance, in [30] - 
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Interestingly enough, at finite temperature, the imaginary part of the effective 
Lagrangian turns out to be independent of the temperature and density and coincides 
with the massless planar limit of the celebrated Schwinger formula [35]. The 
corresponding probability of vacuum decay, for a realistic graphene sample, can be 
very high for laboratory fields, at variance with the common case of massive charged 
electron positron pairs. In the case of 3 + 1 dimensions, such independence was shown 
to hold, to the order of one loop, in reference [2]. Our result in 2 + 1 is stronger: 
indeed, our calculations show that the whole effective lagrangian is independent of 
the temperature and of the chemical potential, within the zeta function regularization 
scheme. 

Furthermore, our treatment, based on the zeta function regularization technique, 
allowed us to obtain the (finite) average value for the minimal quantum current density 
and, thus, a minimal quantum conductivity for graphene, which, at variance with 
most theoretical results Q]3 US] agrees, for a particular selection of the phase of 
the determinant in both irreducible representations, with the one exhibited in [5]. 
The corresponding current is topological in origin, as is the behavior of the Hall 
conductivity (for general studies of topological currents in 2+1 QED see, for example, 
(36] [29]). We hope the present paper will be helpful in clarifying the so called "mystery 
of the missing pi" [8] , or its absence thereof [37] . 
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Appendix : Parabolic cylinder functions 



The parabolic cylinder functions, of the special form we are interested in the present 
context, can be defined, e.g., by the integral representation 9.241 1. p. 1092 of ref. [53] 

£>_ iA/2 [±(l + i)£] = _L 2-^/a+i/a e -*A/4 e <e»/2 
V 71 " 

/OO 
- OO 

where A > , (el, argar a / 2 = A/2 for x < , so that 

D* /2 [±(1+<K] = -L 2 <V2+l/2 e -*A/4 e -tf a /2 

i A/2 e -ax a T2te(l-i)C d!B< (2) 



X 

J — oc 

After the change of variable x i — ► — x 

D* /2 [±(l+»)£] = -L 2*V2+l/2 e7 rA/4 e ^/2 

' V 71 " 

a .iA/2 e -2^±2ix(l-i)€ dXy (.3) 
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we eventually come to the conjugation property 

ir a/2 [±(l + i)£] = D iX/2 [±(l-i)a, (.4) 
as naively expected. The following special values appear in our calculations, 

D ±lX/2 (0)=n-^2^T^±^j cosh^, (.5) 

^ ±jA/2 _ 1 (0) = ±-- 1 /2 2± W4-i/2 r (^ ± !A\ s . nh vrA. ( _ 6) 



'±i\/2-i \v) = ' z ■ < i \^r— j sum — 

±^\D ±tX/2 _ 1 (0)\^±smh^, (.7) 

l^ ±l A/2(0)| 2 = cosh^. (.8) 
The parabolic cylinder functions fulfill the recursion formulas 

^-D v (z) = -\zD v (z) + vD u _ l {z), (.9) 
az 2 

-j- D u (z) = - zD v {z) — D v+ x(z) . (.10) 
dz 2 

Consider the combination 

D + = D_ iX/2 [(l+i)(]D iX/2 [(l-i)£] 

+ ^0_a/a-i[(l + iK]^iA/a-;t [(!-»)£]• (-H) 
From the recursion formulae we get 

2^D- i A/a[(l + i)f]^A/2[( 1 - i )f] = 

A(l + i)D_ iA/2 [(H-i)f]I>iA/2- l [(l-*)e] + c.c. 

a ^ D-a/2-i [ (i + i) e] OiA/2-i = 

- A(l + i)£>_ iA /2[(l + i)^]£»a/2-i[(l-i)?] + c.c, 

so that the above combination D + does not depend upon £ and from the conjugation 
property (Q| we can write 

D+= |D_ iV 2(0)| 2 + ^|£>-iA/2-i(0)| 2 = cxp{^A/4}. (.12) 
Let us now consider the further combination 

D- = D_ iX/2 [ (1 + i) £] D lX/2 [-(l-t)C] 

- ^£>-a/2-i[(l+i)C]^ a/2-1 [-(l-i)C]- (-13) 
From the recursion formulae we get 

2 ^ {^a/2 [(*- IK] £>-a/2 [(* + !)£]} = 

A (1 - i) D lX/2 [(i — i)£]D _ a/2 _i [ (i + 1) £] - 

- A (1 + i)D_ iX/2 [ (i + l)(]])a/2-i [(»-!)£] 

A ^ {£>-a/ 2 -i [ (i + 1) £] £»iA/ 2 -i [ (i - i) £]} = 
A (1 - i) D lX/2 [(i — l)£]D _ a/2 _! [ (i + 1) £] - 

- A(l + i) D _ iX / 2 [(« + !)£] -DiA/2-i [(*—!)£], (-14) 
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which leads to the conclusion that also the quantity £>_ is independent of £, and yields 



iA/2l 



i\/2-l\ 



exp{— 7rA/4} . 



(.15) 



The above important properties of the parabolic cylinder functions can be summarized 
in the remarkable formula 



Da 



^ X \D, 



>± = \D iX/2 (0)| 2 ± - |^ a/2 
Consider the second order differential equations 



_i(0)| 2 = exp{±7rA/4}. 



^ + e + X± i)/±(£,A)=0. 



Two pairs of linearly independent solutions for the upper sign equation are 
4 1) (±^,A)= J D_ iA/2 (±^e ri / 



/1 2) (±C,A) 



D 



iX/2-1 



- iri/4 



(.16) 

(.17) 

(.18) 
(.19) 



while two couples of linearly independent solutions for the lower sign equation are 



/ a) (±C,A) 



iA/2-1 



f m {±i,\) = D lX/2 (i^e — 



e 

- 7Tl/4 



/4 



(.20) 
(.21) 



To the aim of verifying linear independence we have to compute the wronskian. Let 
us first calculate derivatives by means of the recursion formulae (O and (j.lOp that 
yield 



— D 



i\/2 

i£D_ iX /2 



±£\/2< 



i/4 



x f±e\/2e m/ ' 



and thereby 



TV 



r sinh — 



ttA 



On the other side we readily find 



W 



(± & A) , f+ z> (± £, A) = T (1 - i) ex P { ^A/4} 



•(2) 



(.22) 
(.23) 



and analogous relationships for the other solutions. 



In order to understand the physical meaning of the solutions of the wave field 
equations, we have to analyze the leading asymptotic behavior of the parabolic cylinder 
functions. Then from eq. 9.246 1. p. 1093 of ref. [26] we have 
D-ix/2 (£%/2eW4) „ (2 £ 2 )-*V4 e ^/8 exp r_ ie 2 /2 } - 
D-a/2-1 (e%/2e"/4) „ 0(r i) 

If instead ( « -A we have either (e"' 4 = | £ | e 57ri/4 or else £e™/ 4 = | £ | e" 37 "/ 4 . 
Now, for arg(£ e 7 ™/ 4 ) = 57ri/4, no reliable asymptotic expansion is available, so that 



U> A>0) 



Planar QED 26 

from eq. 9.246 3., p. 1094 of ref. [55] we obtain the bona fide leading behaviour for 
£ <C — A < : namely, 



D-ix/2-i (£\/2e™/ 4 ) = D_ iA /2-i (| 

2^ /„^iA/4 ___ f ^A , i£ 2 



e -37ri/4 



f2£ 2 V A / 4 exp - — + ^M, (.24) 
r(l + iA/2) v w P \ 8 2 r 



D. iX/2 (\ti\V2e-^) ~ (2£ 2 )- i V4 exp |. 



3ttA ^ 
~~8 2~ 



Of course, the situation becomes exactly time-reversed for the two other linearly 
independent solutions: namely, for ^ > A > we find 

B-iA/a-i (-£V2e-/ 4 ) - D_ iX /2-i (tV2e~ 3 ^ 4 



r(l+»A/2) v s y *\ 8 2 

D-ix/i (-CV2c"/ 4 ) = D_ iA/2 (eV2c- 37ri / 



(2a-V4 exp J_ , (-25) 



;c«-a<o) 



2 

whereas for £ <C — A < we obtain 

£>-iA/2 (|£|^e"/ 4 ) ~ (2£ 2 )- lA /4 e -A/8 exp {„ ze 2 /2 } 
O-iA/2-1 (|£|%/2c-/4) „ 0(r l) 

For a given particle momentum p x = p , we shall associate the stationary asymptotic 
phase 

^ - ^£ 2 (*) = - ^ e ^ 2 +P*-^ ( ' 26) 

to the positive frequency solutions £ 2 (i) which describe a particle, i.e. an electron of 
momentum p , while the stationary asymptotic phase 

px + ^ 2 (t) = P x + \eEt 2 -pt+^ 27 ) 
will describe an antiparticle, i.e. a positron of momentum — p. 
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